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We examine the stiffness of the Heisenberg spin-glass (SG) model at both zero temperature 
(T = 0) and finite temperatures (T ^ 0) in three dimensions. We calculate the excess energies 
at T — which are gained by rotating and reversing all the spins on one surface of the lattice, 
and find that they increase with the lattice size L. We also calculate the excess free-energies 
at T ^ which correspond to these excess energies, and find that they increase with L at 
low temperatures, while they decrease with increasing L at high temperatures. These results 
strongly suggest the occurrence of the SG phase at low temperatures. The SG phase transition 
temperature is estimated to be Tsg ~ 0.19 J from the lattice size dependences of these excess 
free-energies. 
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§1. Introduction 

Spin-glasses have attracted great challenge for computational physics. It has been believed that 
the spin-glass (SG) phase transition occurs in three dimensions (3d) for the Ising model 0!'!) but 
not for the XY and Heisenberg models. HHBB However, numerical studies during the last decade 
have revealed that the SG phase might be more stable than what were previously believed. For 
example, in the isotropic Ruderman-Kittel-Kasuya-Yoshida (RKKY) model, it was suggested that 
the SG phase transition takes place at a finite temperature.!' & It was also suggested that, for 
the Ising model, the SG phase might be realized at low temperatures even in two dimensions 
(2d)E00i) For the XY and Heisenberg SG models, Kawamura and coworkers took notice of 
the chirality described by three neighboring spins, and suggested that, in 3d, a chiral glass phase 
transition occurs at a finite temperature without the conventional SG order. The basis 

of their suggestion is the stiffness of the system at zero temperature.Q'lll) They estimated 9 C > 
and 9 S < from the lattice size dependences of the chiral and the spin domain-wall energies, 
where 9 C and 9 S are the stiffness exponent of the chiralities and that of the spins, respectively. 
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Recently, however, the ambiguity was pointed out for calculating the domain-wall energy. 
In particular, Kosterlitz and Akino (KA) calculated the spin domain- wall energy of the XY 
SG model using a different method and found that 9 S > 0, contrary to Kawamura and coworker's 
suggestion. If KA's result is true, the question arises for the Heisenberg model whether the SG 
phase transition occurs at a finite temperature. 

The purpose of the present paper is to examine the stability of the SG phase in the Heisenberg 
model. We consider the stiffness of the model at zero temperature (T = 0) and finite temperatures 
(T ^ 0) calculating the excess energy and the excess free-energy, respectively, which are gained 
by rotating or reversing all the spins on one surface of the lattice. We find that the excess en- 
ergy and the excess free-energy at low temperatures increase with the lattice size L. These results 
strongly suggest the occurrence of the SG phase at low temperatures. We estimate the transition 
temperature of T$q ~ 0.19 J from the lattice size dependence of the excess free-energy for various 
temperatures. In section 2, after pointing out the ambiguity of the defect energy method conven- 
tionally used to estimate the domain-wall energy, we propose a method for considering the stiffness 
of the system. In section 3, we calculate the excess energy and discuss the stability of the SG phase 
at T = 0. In section 4, we calculate the excess free-energy at T ^ using the Monte Carlo twist 
method© Section 5 is devoted to the conclusions. Preliminary results of the present work were 
reported in Ref. 21. 

§2. Model and Method 

We start with the model described by the Hamiltonian 

ft = " E J a s i ■ s i 

<ij> 

where Si is a Heisenberg spin with \Si\ = 1 and the sum < ij > runs over all nearest neighbor 
pairs. The coupling constant Jy takes +J or — J with the same probability 1/2. 

Let us briefly describe the defect energy method conventionally used so far. The method orig- 
inates from an application of the renormalization-group concept .SB!) That is, one evaluates an 
effective coupling J e g between block spins generated by renormalization. To estimate J^g, one 
considers the domain-wall energy AE, which is usually defined as the difference in the ground state 
energy between two systems, A and B, with the same bond distribution but with different bound- 
ary conditions. For the system A, a periodic boundary condition is applied in every direction, and 
for the system B, an antiperiodic boundary condition is applied in one direction and the periodic 
boundary conditions in the other directions. However, the meaning of the domain-wall energy is 
not clear, because some domain-wall will arise in the system A and some different domain-wall will 
arise in the system B. Then, one might merely examine the difference in the energy between those 
domain- walls. It is also not clear whether the domain- wall energy AE really provides the effective 
coupling J e g between the block spins. 
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Then, apart from the renormalization-group concept, we consider the stability of the spin con- 
figuration of the system itself. The strategy of our study is as follows. We prepare a cubic lattice 
of L x L x (L + 1) with an open boundary condition in one direction (z-direction) and the periodic 
boundary conditions in the other directions. That is, the lattice has two surfaces and f^L+i- 
First we determine the ground state spin configuration of the lattice. Hereafter, the ground state 
spin configurations on fti and are denoted as {Sj^ 1 ^} and {S^ L+1 ^}, respectively. In this spin 
configuration, any distortion (domain-wall) will be removed, because the lattice has free surfaces 
$7i and ^l+i- Then we add a distortion inside the system in the manner that are fixed 

and{Si( i+1 )} are changed under the condition that the relative angles between the spins are fixed. 
The ground state energy of this system is always higher than that of the system before changing 
{S 4 (L+1) }. This excess energy is the net one added inside the system, because the surface energy of 
^L+ii which is given as the sum of the exchange energies between the spins on f^L+ii is conserved. 
We consider the stability of the system on the basis of this excess energy. One might think that 
the fixing of the relative spin directions on Qi and overestimates the stability of the spin 

configuration. We think, however, that this restriction is not serious for discussing the stability, 
because the increase of the excess energy to infinity for L — > oo means nothing but the existence 
of a strong correlation between the spin configurations on Q± and Ql+i- In fact, the same method 
was successfully applied to the Ising SG model in 2dJlHll' 

§3. The Stiffness at T = 

In this section, we consider the stiffness of the system at T = calculating two kinds of excess 
energies. One is the excess energy which is gained by rotating {S^ L+1 ^} and the other is the excess 
energy which is gained by reversing {5i (L+1) }. We think that it is sufficient to examine these two 
excess energies for considering the stiffness, because we can change {S^ L+1 ^} into any direction by 
combining the rotation and the reversal. 

First we consider the stiffness of the system when we rotate {S^ L+1 ^}. The spins {S^ L+1 ^} are 
rotated around a common axis by the same angle <j). Here we note that effective rotating angles of 
the spins on are not the same, because {S^ L+1 ^} have different directions even in the ground 

state. That is, the spin parallel to the axis does not rotate, whereas the spin perpendicular to the 
axis rotates by the angle (p. Then, as a measure of the rotating angle, we introduce the averaged 
rotation angle to defined by 




(3.1) 



where {S i (L+l) {(f))} is the set of the spins on after the rotation. Hereafter, we call the system 

with the surface spin configurations {S^ 1 -*} and {Si^ L+1 ^} the reference system and that with the 
surface spin configurations {S 1 ^ 1 ^} and {Si^ L+1 ^ ((/))} the rotated system. An important point is 
that there is an upper bound of u) which can be reached by this rotation. If {S^ L+1 ^} are randomly 
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distributed, then the maximum value w miB for L — ► oo is given as 




Therefore, it will be sufficient for us to consider the stiffness for uj < o; max , because we consider 
the stiffness for L — > oo. In fact, the maximum value of uj decreases as the number of the spins on 
CIl+i is increased (see Fig. 1). 
The excess energy AE TOt (uj) due to the rotation may be given as 

AE vot {uj) = E rot (uj) - E G , (3.3) 

where Eq and E rot (uj) are the ground state energy of the reference system and the lowest energy 
of the rotated system with uj, respectively. If AE ro t(uj) for uj < w max increases with the lattice size 
L, the ground state spin configuration will be stable. Otherwise, the spin configuration is unstable. 
This concept is, of course, justified in the ferromagnetic Heisenberg model in (^-dimensions, that is, 
AE Iot (uj) = \jui 2 L B , where 9 (= d — 2) is the stiffness exponent. 

In this calculation, it is not easy to obtain the correct value of E ro t(uj), because there are infinite 
rotating axes which give the same angle uj. Here we estimate E ro t(uj) as follows. For each uj, we 
prepare N r (= 200) rotating axes and calculate the minimum energy for each of those axes. The 
energy E rot (uj) is approximated by the lowest value of those minimum energies. Of course, E TO t(uj) 
is not the correct one but an upper bound of it, because N r is finite. In fact, E rot (uj) for N r /2 is 
slightly larger than that presented here (see Fig. 2). However, we think AE ro t(uj) obtained here 
gives qualitatively correct uj- and L- dependences, because the data for N r and those for N r /2 yield 
almost the same values of exponents a and # ro t defined below. 

The lattice sizes studied here are L = 3 ~ 8 and, for each L, the sample averages are taken over 
about 1000 independent bond realizations. In the calculation of Eg and E TOt (uj), we prepare about 
1000 initial spin configurations, which is obtained by annealing the system down to T = 0.1J from 
a higher temperature, and apply the spin-quench algorithm.B& 

In Fig. 1, we present the results of the w-dependence of [A.E ro t(u;)] for various sizes of L, where [■■■] 
means the sample average. We find that for u < uj max [AE rot (uj)] increases with L, implying that 
the ground state is stable against the rotation. To confirm this view, in Fig. 2, we plot [A.E rot (u;)] 
as functions of L for various w in a logdog form. All the data for each uj lie on a straight line with 
almost the same slope 0.8, that is, [A£ , rot (ti;)] a L 9lot with 6> ro t = 0.8 ±0.1 for uj < w max . We also 
examine the w-dependence of the coefficient and find that the data for uj < w max fit the relation 

[AE TOt (uj)} ~Auj a L e ^, (3.4) 
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with a ~ 1.9, 9 TOt ~ 0.8, and A ~ 0.14J. This finding is illustrated by the scaling plot in Fig. 3. 

This result is quite interesting, because it suggests that the SG phase is realized at low tempera- 
tures. A surprising point is that the values of a and 9 m t are close to those of a = 2 and 9 = 1 of the 
ferromagnetic Heisenberg model in 3d, which implies that the response of the spin configuration 
of the SG model against the rotation is analogous to that of the ferromagnetic model. In fact, we 
found that the rotated angles of the spins gradually increase from layer to layer, reminicent of the 
spin- wave excitation.© We note, however, that the coefficient of A ~ 0.14 J is much smaller than 
that of A = 0.5J in the ferromagnetic model. 

Next we consider the stiffness of the system when we reverse {Si (L+1) }. Hereafter, we call the 
system with the surface spin configurations {■S'j'- 1 ^} and {—S^ L+1 ^} the reversed system. The 
excess energy AE rev due to the reversal may also be given as 

AE rcv = -EVev — Eg, (3.5) 

where E xev is the ground state energy of the reversed system. It should be noted here that 
{— S^ L+1 ^} is never realized by any rotation of { < Sj( L+1 )}JlI' Thus A.E rev will give a different 
view on the stability of the system. The sample averages are taken over about 4000 (L = 3 ~ 6) 
and 2000 (L = 7, 8) independent bond realizations. The result of the L-dependence of [Ai£ re v] is 
presented in Fig. 4. We also find that [AE? rev ] increases with L, suggesting that the ground state 
is stable against the reversal. However, it is difficult to estimate the value of the stiffness exponent 
9 rev for the reversal, because the slope of the curve gradually increases with L. Here, we tentatively 
estimate it as 9 rev — 0.4 it 0.1 from the data for L — 6, 7, 8. It should be emphasized that this value 
might not be the one for L — > oo but a lower bound of it. 

These results of the stiffness for both the rotation and the reversal strongly suggest the occurrence 
of the SG phase transition at a finite temperature. It is noted, however, the estimated value of 
#rev ~ 0.4 is considerably smaller than that of # ro t ~ 0.8. There is a speculation that # ro t could not 
be larger than ^rev© Could 9 VCV be equal to 9 rot for L — > oo ? In fact, the result of [AEVev] shown in 
Fig. 4 implies the increase of 9 rev with L. To examine this possibility, we also calculate the defect 
energies [AE^ct] of the system using replica boundary conditions,© where AE^ef = Er — Eq. 
That is, we prepare two replicas with the same bond distribution and two surfaces S\ and Si- The 
surfaces S\ of the two replicas are connected with each other by the ferromagnetic bond of (J, J, J), 
where (J, J, J) is the set of the connected bonds in the three spin components. We calculate the 
ground state energy Eq of the system under the condition that the other surfaces S2 are also 
connected by (J,, J, J). Then we calculate Er under two different boundary conditions for the 
surfaces S^- In the case A, they are connected by the antiferromagnetic bond of (—J, —J, —J), 
which may correspond to the reversal of {S^ L+1 ^} in the present model. In the case B, they are 
connected by an asymmetric bond of (— J, — J, J), which may correspond to the rotation around 
the z-axis by <f> = ir. These results are presented in Fig. 5. We find that, the defect energies in 
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both cases increase with L and their slopes seem to approach as L is increased. However, further 
studies are necessary to examine the possibility of 6 rev = # ro t for L — > oo. Anyway, all the results 
of the stiffness examined here suggest the occurrence of the SG phase at low temperatures. 

§4. The Stiffness at T / 

In this section, we consider the stiffness of the system at T ^ to estimate the SG phase 
transition temperature T$q. We examine two kinds of excess free-energies which correspond to 
the excess energies AE TOt (tv) and AE Tev examined in § 3. First we prepare the reference system 
in which the spin configurations {S'/ 1 -'} and {S^ L+1 ^} on Q\ and f^L+i ar e fixed to those of the 
ground state. Then we calculate the excess free-energies AF TOt (T) and AF Tev (T) which are gained 
by rotating and reversing {S^ L+1 ^}, respectively. Here, we consider AF TO t(T) when {S^ L+1 ^} is 
rotated around a fixed axis (the z-axis) by 7r/2, because it is not easy to find the rotating axis 
which gives the minimum free-energy of the rotated system. We calculate AF rot (T) and AF rev (T) 
using the Monte Carlo twist method proposed by Ueno.0) These excess free-energies are defined 
as 

AF rot , rcv (T) = F rot , rev (T) - Fq (T) , (4.1) 

where Fq{T) is the free-energy of the reference system, and F TOt (T) and F Tev (T) are the free-energies 
of the rotated and reversed systems, respectively. We may rewrite eq.(4.1) as 

A-F rotjI . ev (T) = A£' rotjrev (T) — TAS , roti rev(T), (4.2) 

where AE TO t(T) and AS , rot (T) are the differences in the internal-energy and the entropy between 
the reference system and the rotated system, respectively, and AE rev (T) and AS rev (T) are those 
between the reference system and the reversed system. Here A£ ro t(T) and A.E rev (T) are obtained 
by using the conventional Monte Carlo method. On the other hand, AS TO %(T) and AS Tev (T) are 
obtained by means of the numerical integration of AE TOt (T) and AE Tev (T) from a high temperature, 
that is, 

AS rot ,MT) = [ T ^^'l (4.3) 
J To 

where To is a temperature high enough to satisfy A£ , rot (To) = and AZ? rev (To) = 0. 

The lattice sizes studied here are L = 4 ~ 10 and, for each L, the sample averages are taken over 
about 2000 independent bond realizations. The starting temperatures are Tq/J = 3.0 {L = 4), 
T / J = 2.1 (L = 5), T / J = 1.7 (L = 6), T /J = 1.5 (L = 7), T /J = 1.3 (L = 8), T /J = 1.1 (L = 
9) and Tq/J = 1.0 (L = 10). We use a gradual-cooling method and take 4 ~ 5 x 10 4 MCS for the 
calculation of A£^ ro t(T) and A£' rev (T). 

In Figs. 6 and 7, we show the temperature dependences of [AF rot (T)] and [AF rev (T)], respec- 
tively. The data for T = obtained by using the method in § 3 are also plotted in these figures. 
We find that, at high temperatures [AF ro t(T)] and [Ai ? rev (T)] decrease with increasing L, whereas 
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at low temperatures they increase with L. These results suggest that the SG phase will be realized 
at low temperatures. Usually, one estimates the phase transition temperature from the crossing 
temperature of the free-energies for various lattice sizes L. In the present model, however, we could 
not use this method directly, because the crossing temperature Tl for the lattice sizes L and L + 1 
shifts systematically to the low temperature side with increasing L. Then, we tentatively assume 
that Tl's for [AT rot (T)] and [AT rev (T)] decrease linearly with 1/L, and plot them as functions 
of 1/L, which are shown in Fig. 8. In fact, Tl for each of the cases seems to lie on a straight 
line. Using the least-squares method, we estimate T^/J = 0.192 ± 0.015 for the rotation and 
Too/ J = 0.188 ± 0.015 for the reversal. A remarkable point is that these values of Tl for L — > oo 
are almost the same. Therefore we may conclude that, if the SG phase transition really occurs, the 
transition temperature should be Tsg ~ 0.19 J. 

§5. Conclusions 

We have examined the stiffness of the Heisenberg SG model in 3d at both T = and T / 0, 
using a method which concerns with the stability of the spin configuration of the system itself. We 
calculated the excess energies [Ai? rot (w)] and [AS rev ] at T = which are gained by rotating all the 
spins on one surface of the lattice and reversing them, respectively, and found that both [A£' rot (a;)] 
and [AE'rev] increase with the lattice size L. This result strongly suggests the occurrence of the SG 
phase transition at Tsg 0- Then, we calculated the excess free-energies [AT rot (T)] and [AF rev (T)] 
at T ^ to estimate the transition temperature. We found that each of these free-energies increases 
with L at low temperatures, whereas it decreases with increasing L at high temperatures. From 
these results, we estimated the SG phase transition temperature as Tsg ~ 0.19 J. Recently we also 
performed the Monte Carlo simulation and found an evidence of the SG phase transition at almost 
the same temperature of Tsg ~ O.18J.0 Hence we believe that the SG phase is really realized in 
this model and the transition temperature to this phase is T§q/J = 0.19 ± 0.02. However, further 
studies are required to confirm our view. We hope that the present work stimulates other works 
on vector spin glasses. 
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Fig. 1. The excess energies [A£ rot (w)] of the Heisenberg SG model in 3d for various lattice sizes of L as functions 
of the rotating angle w. The arrow indicates the position of w max - 

Fig. 2. The lattice size dependences of the excess energies [Ai? rot (w)] of the Heisenberg SG model in 3d for various 
rotating angles uj. Cross and circle represent the data for 100 and 200 rotating axes, respectively. 

Fig. 3. Scaling plots of the excess energies [A£ ro t(w)] of the Heisenberg SG model in 3d. Here Io t = 0.8 and 
a = 1.9. We show the results only in the uj- range which is relevant in the model (the arrow indicates the position 

Of Wmax)- 

Fig. 4. The lattice size dependence of the excess energy [AE lav ] of the Heisenberg SG model in 3d. 

Fig. 5. The lattice size dependences of the defect energies of the Heisenberg SG model in 3d which are calculated 
by using two replica boundary conditions explained in the text. 

Fig. 6. Temperature dependences of the excess free-energies [AF Iot (T)] of the Heisenberg SG model in 3d. Note 
that [AF rot (T)] at T = is the excess energy [Ai? rot ] which is gained by rotating around the z-axis by tv/2, and it 
is calculated by the different method in § 3. 

Fig. 7. Temperature dependences of the excess free-energies [AF rBV (T)] of the Heisenberg SG model in 3d. Note 
that [AF rev (T)] at T = is the excess energy [Ai5 rcv ] calculated by the different method in § 3. 

Fig. 8. The 1/L dependences of the crossing temperature Tl of the excess free-energies for the lattice sizes of L 
and L + 1 . 
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